On the quantum discord of general X states

Quantum discord Q is a function of density matrix elements. The domain of such a function in the case of two-qubit system with X density matrix may consist of three subdomains at most: two ones where the quantum discord is expressed in closed analytical forms (Q π/2 and Q 0) and an intermediate subdomain for which, to extract the quantum discord Q θ , it is required to solve in general numerically a one-dimensional minimization problem to find the optimal measurement angle θ ∈ (0, π/2). Hence the quantum discord is given by a piecewise-analytic-numerical formula Q = min{Q π/2 , Q θ , Q 0 }. Equations for determining the boundaries between these subdomains are obtained. The boundaries consist of bifurcation points. The Q θ subdomains are discovered in the generalized Horodecki states, in the dynamical phase flip channel model, in the anisotropic spin systems at thermal equilibrium, in the heteronuclear dimers in an external magnetic field. We found that transitions between Q θ subdomain and Q π/2 and Q 0 ones occur suddenly but continuously and smoothly, i.e., nonanalyticity is hidden and can be observed in higher derivatives of discord function.